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ABSTRACT 

We present constraints on the non-linear coupUng parameter /„; with the Wilkinson Microwave 
Anisotropy Probe (WMAP) data. We use an updated method based on the spherical Mexican hat 
wavelet (SMHW) which provides improved constraints on the /„; parameter. This paper is a continu- 
ation of a previous work by Curto et al. where several third order statistics based on the SMHW were 
considered. In this paper, we use all the possible third order statistics computed from the wavelet 
coefficient maps evaluated at 12 angular scales. The scales are logarithmically distributed from 6.9 
arcmin to 500 arcmin. Our analysis indicates that /„; is constrained to —18 < /„; < -1-80 at 95% 
confidence level (CL) for the combined V-l-W WMAP map. This value has been corrected by the 
presence of undetected point sources, which adds a positive contribution of A/„i = 6 ± 5. Our result 
excludes at ~99% CL the best-fitting value /„; =87 reported by Yadav & Wandelt. We have also 
constrained /„; for the Q, V and W frequency bands separately, finding compatibility with zero at 95 
% CL for the Q and V bands but not for the W band. We have performed some further tests to un- 
derstand the cause of this deviation which indicate that systematics associated to the W radiometers 
could be responsible for this result. Finally we have performed a Galactic North-South analysis for 
fni- We have not found any asymmetry, i.e. the best-fitting /„; for the northern pixels is compatible 
with the best-fitting /„; for the southern pixels. 

Subject headings: methods: data analysis - cosmic microwave background 



1. INTRODUCTION 

The cosmic microwave background (CMB) offers a pic- 
ture of the early Universe when it was only 400,000 
years old. The CMB photons last scattered off electrons 
at that time and since then they traveled free through 
the space. The primordial perturbations set up dur- 
ing infiation are imprinted in both radiation and mat- 
ter distribution. The CMB temperature anisotropics, 
related with the primordial perturbations, can be used 
to test some assumptions of the so called standard 
model. In particular, we can test the pred iction of the 
standard, single field, slow roll in f iation (iGuthI 119811 : 
lAlbrecht fc SteinhardU [1981 [Linda [1981 . 1198.'^ which 
states that the anisotropics are compatible with a nearly 
Gaussian random field. Other non-standard models of 
inflation predict detecta ble levels of non-Ga ussianity in 
the anisotropics (see e.g. lBartolo et al.ll2004D . 

Primordial non-Gaussianity of the local form is char- 
acterised by the non-line ar coupling paramet er /„,; 
jSalopck fc Bond|[l9 90: Ga ngui et al.|[l994: Verde et all 
[20001: IKomatsu fc Sperge]|l200l[k 



the angular bispectrum and wavelets on the Cosmic 
Back gr ound Explorer (C OBE) data (jKomatsu et al.l 
l2002t iCavon et al.l 120031). the angular bispectrum 
on MAXIMA data (jSantos et all 120031) . the angu- 
lar bi spectrum on WMAP data (IKomatsu et al.l 



2006 



20031: iCreminelli et al.l , 

2007t lYaday k Wandelt [200l IKomatsu et al. 200L, 
Smith et al.ll20p^. different kind of wave l et analyses on 
WMA P data (IMukheriee fc Wand [20f)l ICabella et all 
120051: iCurto et al.l 120091 ), the Minkowsk i func tionals 
BOOMERanG data pe Troia et al.l I2007D. the 



on 



<i>(x)=<i>i(x) + /„,{$i(x)-(<j>i(x))} 



(1) 



where $(x) is the primordial gravitational potential and 
$l(x) is a linear random field which is Gaussian dis- 
tributed and has zero mean. 

There are many studies based on different statis- 
tical tools to constrain the local /„; parameter from 
the CMB anisotropics, using the data of different 
experiments. We can mention the analyses using 
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Minkowski functionals on Archeops data (jCurto et al.l 

fD7, 2008 ), the Mink o wski functiona l s on WMAP data 
omatsu et al.l 120031: iSoergel et al.l 120071: iGott et all 
120071: iHikage et al.l 120081 : IKomatsu et al.l 120091 ) among 
others. We can also mention new promising techniques 
as for example one based on th e n-point probability 
density distributio n (iVielva fc S anz 2009), and other 
based on needlets (jPietrobon et al. 2009: Rudjord et al.l 
l2009f ). Other works are based on the use of the large 
scale structure to constrain (jSlosar et al.ll2008[ ). 

This work is a con t inuati on of the wavelet-based anal- 
ysis by ICurto et all ([2009^1 of the WMAP^ data. We 
use high resolution WMAP data maps, and compute 
the wavelet coefficients for 12 angular scales logarithmi- 
cally spaced from 6.9 arcmin to 500 arcmin. With these 
wavelet coefficients we compute all the possible third or- 
der moments involving these scales. 

The article is organised as follows. Section [2] presents 
the estimators used to test Gaussianity and to constrain 
fni, the data maps and the simulations. Section [3| sum- 
marises the main results of this work and the conclusions 

^ http://map.gsfc.nasa.gov/ 
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are in Section 01 

2. METHODOLOGY 

This analysis is based on t he spherical Mexican hat 
wavele t (SMHW) as defined in lMartmez-Gonzalez et ahl 
(|2002f ). For references about the use of the SMHW to 
tes t Gaussianity in the CMB see for example the review 
by iMartinez-Gonzaled (|2008f ). We compute the wavelet 
coefficient maps at several scales Ri logarithmically sep- 
arated {Ri+i/Ri constant). The considered scales are'^: 
Ri = 6.9', i?2 10.6', i?3 = 16.3', i?4 = 24.9', 
i?5 = 38.3', Rq = 58.7', R^ = 90.1', i?8 = 138.3', 
i?9 = 212.3', Rw = 325.8', Ru = 500'. We also in- 
clude the unconvolve d map, which will be represented 
by the scale Rq as in lCurto et al.1 (|2009f) . For each pos- 
sible combination of three scales Ri, Rj, and Rk (where 
the indices i, j, and k can be repeated) we define a third 
order statistic 

Qijk = -r > (,2) 

JVij,fc cTiajak 

p—O 

where Npix is the total number of pixels of the map, 
Ni j^k is the number of pixels available after combining 
the extended masks corresponding to the three scales 
Ri, Rj and Rk, Wp^i — Wp{Ri) is the wavelet coeffi- 
cient in the pixel p evaluated at the scale Ri, and at 
is the dispersion of Wp^i. Each map Wp^i is masked out 
with the corresponding extended mask at the scale Ri 
as in iCurto et al.l |2009) . For a set of n scales we have 
nstat = (n -I- 3 — l)!/[3!(n — 1)!] third order statistics 
such as the one defined in Eq. We have tested with 
simulations that these statistics have a Gaussian-like dis- 
tribution. We can construct a vector q of dimension Ugtat 

q = [90,0,0; 9o,o,i; 90,0,11; 90,1,1; •■•; 911,11,11]- (3) 

With this vector we can perform two different analyses 
using a statistic: one to test Gaussianity and a second 
one t o constrain /„; (Eqs. (7) and (8) of ICurto et al.l 
[2OOI . 

We use the 5-yr WMAP foreground reduced data, 
available in the Legacy Archive for Microwave Back- 
ground Data Analysis (LAMBDA) web site'*. We com- 
bine the maps of different radiometers us ing the inverse 
of the noise variance as an optimal weight (jBennett et al.l 
|2003). In particular we analyse the V-l-W, Q, V and W 
combined maps at a resolution of 6.9 ar cmin, correspond- 
ing to a HEALPix (|G6rski et al.ll2005ll N^tde = 512. We 
use the KQ75 mask and also a set of extended masks for 
the wavelet coefficien t maps. We u se the same masks as 
the ones described in lCurto et al.l (|2009| ) for a threshold 
of 0.01. This corresponds to an available fraction of the 
sky from 71.2% for the i?i scale to 31.4% for the Ru 
scale. Notice that larger scales have the restriction of a 
lower available area, which means a lower sensitivity to 

fnl • 

Finally we analyse the data and compare them with 
Gaussian and non-Gaussian simulations. The Gaussian 

^ The level of discretisation of the wavelet space is a balance 
between the minimum number of scales needed to extract the non- 
Gaussian signal from the data and acceptable computational re- 
quirements. 

* http://lambda.gsfc.nasa.gov 



TABLE 1 

CONSTRUCTED FROM THE 364 STATISTICS FOR V+W, Q, V, 
AND W DATA MAPS. We ALSO PRESENT THE MEAN AND THE 
DISPERSION OF THE CORRESPONDING TO 1,000 GAUSSIAN 
SIMULATIONS, AND THE CUMULATIVE PROBABILITY FOR THE X^ OF 
THE DATA OBTAINED FROM THE SIMULATIONS. 



Map 




DOE 


ix') 


(J 


P{X' < xLta) 


V-l-W 


349 


364 


379 


49.2 


0.29 


Q 


384 


364 


378 


48.7 


0.63 


V 


348 


364 


377 


47.8 


0.27 


w 


354 


364 


376 


47.9 


0.35 



TABLE 2 

Best-fitting /„; values obtained from V+W, Q, V and W 

COMBINED maps. We ALSO PRESENT THE MEAN, DISPERSION AND 
SOME PERCENTILES OF THE DISTRIBUTION OF THE BEST FIT /„; 
VALUES OBTAINED FROM GAUSSIAN SIMULATIONS. 



Map 


best /„; 


{fnl) 




^0.160 


^0.840 


^0.025 


^0.975 


V-l-W 


39 


-1 


25 


-26 


24 


-51 


47 


Q 


11 





33 


-31 


34 


-63 


66 


V 


23 





30 


-28 


30 


-55 


59 


w 


65 


-4 


30 


-33 


26 


-59 


58 



simulations are performed using the best fit power spec- 
trum d for WMAP provided by LAMBDA and the 
instrumental white noise of each WMAP radiometer. 
The non-Gaussian simulations with the /„; contribu- 
tion are co mputed following the algorithms described in 
iLiguori et al. (2003, 200"^) and transformed into WMAP 
maps with the instrumental noise included. We also es- 
timate the unresolved point source contribution to /„; 
for the V-j-W case by analysing point source simulations. 
These simulations have been generated as in lCurto et al.l 
([2009) foll owing t he source number counts dN/dS given 
bv ide Zotti et all (|2005D . 

3. RESULTS 

In this Section we present the Gaussianity analysis of 
the WMAP data for the combined V+W, Q, V, and W 
data maps. We also constrain /„; for these maps through 
a test. We estimate the contribution of point sources 
to the V+W map. Finally we constrain /„; for northern 
and southern pixels separately. 

3.1. Analysis of WMAP data 

We evaluate the wavelet coefficients at the 12 consid- 
ered scales. With these coefficients we compute the third 
order estimators defined in Eq. For 12 scales, we 

have 364 possible third order statistics. We compute 
these statistics for the data maps and for Gaussian sim- 
ulations. The covariance matrix used in the statistics 
is constructed from 10,000 Gaussian simulations. For the 
considered cases, V-l-W, Q, V and W, we have that the 
data are inside the 2a error bars, i.e., the data are com- 
patible with Gaussian simulations. We compute a x^ 
statistic by comparing the data wit h the expect e d valu e 
of Gaussian simulations following ICurto et al.l (|2009D . 
We as well compute the of an additional set of 1,000 
Gaussian simulations. The statistic of the data is 
compatible with the of the Gaussian simulations. This 
is presented in Table [1] 

We have studied the properties of the covariance ma- 
trix for the V-l-W map (see left panel of Fig. [T]) which is 
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Fig. 1. — Left panel: the covariance matrix elements of the third order statistics for 
Cijk,rst/ \/Cijk,ijkCrst,rst > 0.95. Right panel; the covariance matrix eigenvalues for the V+W map. 
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Fig. 2. — The expected values of the normalised third order statis- 
tics for V+W simulations with different values of where 
r = k} is ordered following Eq. (|3]l. 
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Fig. 3. — The normalised third order statistics qr for the com- 
bined V+W data map. The solid line corresponds to the expected 
values for the best-fitting /„; model (/„; = 39). We also plot the 
f(T error bars in green. 

derived from 10,000 Gaussian simulations. We may won- 
der if that number of simulations is enough to achieve 
the convergence. The condition number, defined here as 
the ratio between the maximum and minimum eigenval- 
ues (see right panel of Fig. [T]) is Cond{C) ^ 10 . An 
upper limit for the relative errors in the inverse covari- 
ance matrix is CondiC) X Crnach i whcrc f:rnach is the Com- 
puter error, €mach ^ 10^^^ for double precision. Thus 
the relative error for the inverse matrix coefficients is 
Cond{C) X Cmach ~ 10~^ which is almost negligible. We 



have also computed the covariance matrix with two in- 
dependent sets of 5,000 simulations and constrained the 
best-fitting /„; of the data with these two matrices. The 
best-fitting /„; value of the data is almost the same for 
both matrices and it agrees with the result obtained with 
10,000 simulations suggesting that the convergence is al- 
most reached with 5,000 simulations. 

We impose constraints on /„; through a test. We 
calculate the expected values of the estimators for dif- 
ferent fni cases using a set of 300 non-Gaussian simula- 
tions. In Fig. [2] we plot the expected values of the 364 
statistics for several /„; cases for the V+W map. We 
have checked that the estimator is unbiased. We have 
evaluated the expected values of the estimator using 200 
non-Gaussian simulations of the V+W map and analysed 
the remaining 100 independent non-Gaussian simulations 
with /„; — 50. The result is an average best-fitting /„; 
value of 50.6. 

In Fig. [3] we plot the values of the statistics for the 
data map and compare them with the expected values 
for the best-fitting /„; model. Then we perform a 
analysis to find the best-fitting value of /„; for each map. 
We also analyse Gaussian simulations in order to obtain 
the frequentist error bars. Table [2] lists the best-fitting 
/„; values for the V+W, Q, V and W combined maps 
and the main properties of the histograms of the best- 
fitting /„; obtained from Gaussian simulations. In the 
left panel of Fig. d] we plot the x^ifni) versus /„; for 
the V+W map, and in the right panel of Fig. |4]we plot 
the histogram of the best- fitting /„; of 1,000 Gaussian 
simulations. We have -12 < /„; < +86 for V+W, -52 < 
/„/ < +77 for Q, -32 < f^i < +82 for V and +6 < /„, < 
+ 123 for W (aU at 95% CL). Note that fni increases as 
the frequency grows from Q to V and W bands. This 
suggests the possible presence of foregrounds residuals 
as they are more important at low frequencies and they 
add a negative contribution to fn i for the bispectrum 
estimator (jYadav fc Wandeltl 12008( 1. To further check if 
this is also the case for the wavelet estimator, we have 
also studied the K and Ka bands, where the foreground 
signal has an important contribution. We have obtained 
fK ^ _497 ± 42 and /,^" = -18 ± 37. Taking this into 
account and the values of /„; for the clean and raw Q, 
V, and W maps (see Tables [2] and [3]) we can see that 
our estimator is sensitive to the presence of foregrounds 
(biasing the result towards lower values). The sensitivity 
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Fig. 4. — The x^ifnl) statistics versus /„; for the V+W data i 
Gaussian V+W simulations. The dispersion is cr{fnl) = 25. 

is even more significant for the Q, V, and W bands u sing 
the bispectrum (see Table 6 in lKomatsu et al.ll2009f ). 

It is interesting to point out that the wavelet-based 
method has an intermediate dimension (ustat — 364 for 
12 scales) when compared with the bispectrum combi- 
nations 

i^tax where 2 < /? < 3 and i^ax ~ 10^) and 
the Minkowski functionals (usually several tens of com- 
binations). The number of wavelet cubic combinations 
increases significantly with the number of scales. 

We also estimate the contribution of the point source s 
for the V+W combined map as in ICurto et all (pOffflh . 
We add the point source simulations to the CMB plus 
noise simulations. For each one of them we compute its 
best-fitting /„; and compare it with the obtained for the 
same case without including the point source simulation. 
The difference returns an estimate of the contamination 
on fni due to point sources. For the V+W map we have 
^fni = 6±5. Therefore our estimate taking into account 
the point sources is — 18 < /„; < +80 at 95% CL for the 
V+W map. Compar ing with the best-fitting v alue for the 
V+W map given bv lYadav fc WandelH (l2008h . fni = 87, 
our analysis excludes that value at ^99% CL. H owever, it 
is important to point out that lYadav fc Wand elt (2008) 
used different choices for the data maps in the analysis 
(they used WMAP 3-yr data whereas we use WMAP 5- 
yr data, different weighting for the channels, etc.), which 
could also contribute to the found discrepancy. In par- 
ticular, a simple average of the channels (instead of a 
noise weighted combination) enhances the signal at high 
multipoles at the cost of having noisier maps. 

3.2. Constraints on fni with the W band 

The W map best-fitting /„; value is only compatible 
with zero at 99% CL. This result is in apparent dis- 
crepancy with the values obtained for the V and V+W 
maps, which are compatible with zero at 95%CL. The 
point sources add a low contribution to the W map, 
^fni = 1 ± 2, and therefore they do not explain its best- 
fitting value of /„; =65. We may wonder if that value 
can be obtained by a statistical fluctuation. Considering 
simulations with different models (/„; = 0, /„; — 40 and 
fni = 70) we have confirmed that the best-fitting fni for 
W is compatible with the values obtained for the V and 
V+W maps. 

To understand the relatively large /„; value found in 
the W map we have performed some additional tests. 
First of all, we have checked if this deviation could be 
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due to the presence of residual foregrounds by studying 
the V-W map for the clean and the raw (before tem- 
plate subtraction) maps as well as the raw maps for the 
Q,V,W and V+W cases (see Table^. We find that there 
is a very significant (positive) deviation in the clean V- 
W map. Notice that the V-W map has CMB residuals 
due to the different resolutions of the V and W bands, 
which have also been taken into account in the V-W 
simulations. In any case, the signal due to the CMB 
is very small and we do not expect it to affect the results 
given in Table [3] for the different maps. Since this com- 
bination contains mainly residual foregrounds and noise, 
both could be responsible for the deviation. Interest- 
ingly, when we repeat the test for the raw V-W map, 
where foreground contamination should be more impor- 
tant, the best fni value becomes compatible with the 
simulations. This indicates again that foreground emis- 
sion tends to bias the estimated /„; towards lower values. 
This is also observed for the best fni value estimated for 
the raw Q, V, W and V+W maps, which is systematically 
lower than the one obtained for the clean maps. There- 
fore, for the case of the raw V-W map, some effect from 
systematics may be cancelled by foreground residuals. 

If foregrounds are not responsible for the deviation 
found in the W and V-W maps, we may wonder if it 
is due to systematics present in the W radiometers. To 
test this possibility, we have studied two different com- 
binations of the four W radiometers, where CMB and 
foregrounds are basically cancelled. One of these combi- 
nations is consistent with Gaussian simulations but the 
second one shows again a deviation at the 95 per cent 
CL, indicating the possible presence of some spurious 
signal in the noise of one or several of the W radiome- 
ters. In order to localise further the origin of this signal, 
we have also studied each W radiometer separately, find- 
ing a deviation at the level of 98 per cent for the W2 
radiometer, with a best /„; value of 91, while the rest of 
the radiometers are consistent with the zero value at the 
95 per cent CL. Finally we have also considered the dif- 
ference between the two V radiometers, which is found to 
be compatible with Gaussianity. This further indicates 
that, if a systematic is responsible of the V-W map, this 
would be present in the W frequency channel (see Table 

El). 

All these tests suggest that the relatively large /„; 
value obtained for W may come from systematics present 
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TABLE 3 

Best-fitting /„; values for different WMAP radiometers 

AND combinations OF THEM FOR RAW AND CLEAN DATA. We ALSO 
PRESENT THE MEAN AND THE DISPERSION OF THE BEST-FITTING 
/„; VALUES OBTAINED FROM GAUSSIAN SIMULATIONS. 
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foreground 
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^ \Jnl 1 
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raw 


34 


-1 


25 




Q 




raw 


-3 





33 
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raw 


16 





30 
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raw 


60 





30 




V-W 




raw 


-0.02 


0.02 


0.29 




v-w 




clean 


1.02 


0.02 


0.29 




Wi 




clean 


39 
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41 




W2 




clean 


91 


-3 


45 




W3 




clean 


23 


2 


47 








clean 


59 





44 


Wi 


+ W2-W3 


- W4, 


clean 


-0.52 


0.00 


0.34 


Wi 


-W2 + W3 


- W4 


clean 


0.85 


0.00 


0.35 




V1-V2 




clean 


-0.01 


0.01 


0.34 



in the W radiometers. In any case, a more exhaustive 
study is necessary in order to estabhsh the origin of this 
deviation. 

3-3. fni for the North and South hemispheres 

The locahzation property of the wavelets aUows a lo- 
cal analysis of the /„; parameter. In particular, we test 
the Gaussianity and estimate the best-fitting /„; value 
of the V-l-W combined map using only northern (Galac- 
tic latitude & > 0) and southern pixels {h < 0) in Eq. 
[21 For both cases the third order statistics obtained 
for the data are compatible with Gaussian simulations 
(inside the 2(t error bars) . In Table [5] we list the best- 
fitting /„; values for the northern and southern hemi- 
spheres and the main properties of the distribution of 
the best-fitting /„; obtained from 1,000 Gaussian simu- 
lations. We estimated the contribution of the unresolved 
point sources as in the previous subsection, and the val- 
ues are Afni = 7 ± 7 for the North and A/„; = 5 ± 7 
for the South. Taking into account this, the results are 
-32 < /„; < 113 for the North and -50 < /„; < 99 for 
the South at 95% CL. These values are compatible with 
zero at 95% CL. We also study the compatibility of the 
best-fitting value for the North and South hemispheres 
between them. We compute the difference of the best- 



fitting fni value Afnl 



AN) 
Jnl 



f^i'' for the North and 



South hemispheres for the set of 1,000 V-l-W Gaussian 

simulations. The difference is A/^^"*"^ — 11 for the data, 
and for the simulations is A/„i = 3 ±55. The cumulative 

probabihty is P(A/„, < A/^^f"*"') = 0.57 and therefore 
the difference for the data is compatible with the results 
obtained from simulations. This means that we do not 
find any asymmetry in the North- South /„; value. 

4. CONCLUSIONS 

We have tested the Gaussianity and constrained the /„; 
parameter with the 5-yr WMAP data. We use an opti- 
mal wavelet-based test. We have considered the V+W, 
Q, V and W combined maps at high resolution. We 
have used a set of 300 realistic non-Gaussian simulations 
and thousands of Gaussian simulations for the analysis. 
We have computed the wavelet coefficient maps at scales 
from 6.9 arcmin to 500 arcmin and computed all the 



TABLE 4 

Best-fitting /„; values obtained for the northern and 

SOUTHERN hemispheres. We ALSO PRESENT THE MEAN, 
dispersion and SOME PERCENTILES OF THE DISTRIBUTION OF THE 
BEST FIT fni VALUES OBTAINED FROM GAUSSIAN SIMULATIONS. 



Region 


best /„; 


ifnl) 




-'f 0.160 


^0.840 


^0.025 


^0.975 


North 


46 


2 


37 


-35 


40 


-71 


74 


South 


35 


-1 


38 


-39 


37 


-80 


69 



possible third order moments (Eq. [J) using appropriate 
extended masks. 

The data are compatible with Gaussian simulations for 
the considered combined maps (see Table [Ij . We have 
imposed constraints on the non-linear coupling parame- 
ter /„; by using non-Gaussian simulations with The 
results show that /„; increases when we go from the Q 
to the V and W combined maps. Thi s frequency depen- 
dence a.lso appears in t he re s ults bv lYadav fc Wandeltl 
(|2008[ ): [Kmnatsu et all (|2009l ): iCurto et all (|2009D . The 
results are compatible with zero at 95% CL for the V-l-W, 
Q, and V combined maps, but not for the W map (which 
is compatible at 99% CL). This value cannot be explained 
by unresolved point sources since their contribution is 
A/„i = 1 ± 2 for the W map. We have estimated the 
probability of having those values with simulations and 
the results do not show incompatibility among different 
channels. We have also seen that the relatively large /„; 
value obtained for the W band may come from system- 
atics in one or several radiometers of this band. 

We have also estimated the contribution of unresolved 
point sources to /„; for the V+W ma p using a reahstic 
model given by Ide Zotti et al.l ()2005D . The results are 
A/„; =6 + 5. Taking into account this value, our best 
estimate for /„; is -18 < /„; < +80 at 95% CL. The use 
of new scales and all the possible third order moments 
has returned better constrai ns to /„; and lower error bars 
compared with the results bv lCurto et al.l (2009) and pre- 
vious works. Our best estimate is compatible with the 



vious works. Uur best estimate i s comp atible witn tne 
values obtained bv lKomatsu et al.l (|2009|) and it excludes 
the best-fitting /„; value obtained bv lYadav fc Wandeltl 
(2008) at the -99% CL. 

Finally we have constrained /„; for the North and 
South hemispheres and the results give two best-fitting 
values that are compatible with zero at 95% CL and also 
are compatible between them. Therefore we do not find 
any North-South asymmetry for this parameter. 
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